Towards a common thread in Complexity: an accuracy-based 

■^- . approach 

O 

psj . Pragya Shukla 

Department of Physics, Indian Institute of Technology, Kharagpur, India. 

in ' (Dated: July 25, 2008) 

(N' 

Abstract 

O ' The complexity of a system, in general, makes it difficult to determine some or almost all matrix 

elements of its operators. The lack of accuracy acts as a source of randomness for the matrix 

C^ . elements which are also subjected to an external potential due to existing system conditions. The 

fluctuation of accuracy due to varying system-conditions leads to a diffusion of the matrix elements. 
We show that, for single well potentials, the diffusion can be described by a common mathematical 



a 



'^ . formulation where system information enters through a single parameter. This further leads to a 

O . characterization of physical properties by an infinite range of single parametric universality classes. 
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INTRODUCTION 



For systems that can be described mathematically, physical information can be derived, 
in principle, from a detailed knowledge of the operators that govern their evolution. Physical 
systems can however be complex in nature and it is not always possible to determine the 
operator exactly or, even if they are known, to solve the equations they determine. This 
paper aims to model the statistical behavior of those complex systems where a matrix 
representation of the operators is meaningful. 

The complexity may appear in various forms, for example, as noise due to many body 
interactions or external disorder potential, as chaos due to scattering of a particle from 
boundaries (e.g. clean quantum dots), as coherence patterns emerging out of randomness 
(see, for example, l[] for various definitions of complexity). For example, consider the 
Hamiltonian of a many body system . If the local interactions are complicated in a specific 
part of the system, the evaluation of the corresponding matrix elements becomes technically 
difficult. These elements can then be determined only within a certain degree of accuracy and 
can best be described by a probability density. However the system may also contain parts 
where interactions are simple and the related matrix elements can exactly be calculated. The 
operator then turns out to be a matrix with both random and non-random elements; we 
refer such a matrix as a generalized random matrix. Simlar matrices would also appear for 
systems containing a combination of chaotic as well as ordered components. The properties 
of such system can then be modeled by an ensemble of generalized random matrices. 

In recent years, due to increasing degree of complexity in systems of industrial and tech- 
nological interests, the mathematical models such as random matrix ensembles have become 
„ecessa,v. I. fac, a pa.«cu.a. c.ass of *ese ensembles. Unowu as stat,o„ary e„se.bles Q. 
have been successfully applied for modeling of the operators for a wide range of complex 
systems e.g. nuclei, atoms, molecules, disordered and chaotic systems, quantum chromo- 
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dynamics, elastomechanics, electrodynamics (see reviews [2|, 13|, 14, 15|, 16|, 17|, 18|, 19|, 110|, lUH and 
references therein for details), mathematical areas such as Riemann zeta function, enu- 
meration problems in geometry and fluctuations in random permutations [l2|, biological 



systems pjj], stock markets |l4], atmospheric sciences ISj, complex networks [la] 

also [17i]). The stationary random matrix ensembles are basis-invariant ensembles, charac- 



etc (see 



terized by similar and independent distribution of almost all elements [2|. This restricts 
their applicability only to the generators with wave-functions extended in the entire system 
or with a coherent scattering of waves. However the matrix elements distribution can sig- 
nificantly be affected by various system conditions e.g missing interactions among some of 
the sub-systems, a variation in their degree or nature, symmetry and boundary conditions. 



dimensionality, disorder etc 22]. These condition may result in different strengths of the 
elements, correlations between them and localized waves; the corresponding ensembles are 
then basis-non- invariant. 

The presence of local interactions and wave-localization phenomena is quite generic to a 
wide range of complex systems. The statistical analysis of their physical properties requires 
therefore a search for new mathematical tools. The present study is an attempt in this 
direction. The basic idea here is to take into account the inaccuracy in the matrix repre- 
sentation of an operator of a complex system. The fluctuation of accuracy with changing 
system conditions results in a change of distribution parameters of various matrix elements. 
This leads to a seemingly multi-parametric diffusion of the ensemble density. However as 
shown here, the diffusion is essentially governed by a single parameter only. The informa- 
tion can then be utilized to study the fluctuations of the physical properties due to varying 
system conditions and express the results in a common mathematical form for a wide range 
of complex systems. 

The statistical behavior of complex systems and the possibility of a common mathemati- 



cal formulation was recently studied by considering their maximum entropy models 19|, l20] 



The latter are based on the formulation of the ensemble density by maximizing the infor- 



mation entropy under constraints imposed on the system [18|]. The ensemble density is then 
utilized to extract the distribution of eigenvalues and eigenfunctions and desired physical 
information. The maximum entropy approach indicated the possibility of a classification of 
the complex systems into various universality classes (based on the behavior of their statis- 



tical measures and characterized by the complexity parameter 19|, l20|). These results were 



also verified numerically for certain cases 2l|, |22, l26|]. However the complexity parameter 



formulation within this approach gives rise to some queries which required a more intuitive 
physical reasoning for their resolution. This motivates us to consider the accuracy based 
approach which not only resolves the queries but also helps in generalization of the single 
parametric formulation to a wider range of complex systems (those with system conditions 
subjecting matrix elements to a potential with a single minima only). 

The paper is organized as follows. Section II describes the diffusive dynamics for the 
matrix elements of a Hermitian operator subjected to an external potential of type e"^*-"^-* 
as well as random noise originated in complexity of the system. (We have considered here 
the real-symmetric case only however the results are valid for complex Hermitian and real- 
quaternion cases too). The comparison of this approach with maximum entropy approach 
is discussed in section HI. This is followed by section IV describing the derivation of the 
statistical measures of the eigenvalues and eigenfunctions using standard perturbation the- 
ory. Note, for generic potentials, the derivation of the measures by a direct integration of 



the evolution equation (the method used for Gaussian cases in 21j) is technically difficult. 
We conclude in section V with a summary of our main results. 

II. ACCURACY DRIVEN DIFFUSION OF MATRIX ELEMENTS 

Consider, as an example, a Hermitian operator H oi a. complex system with time- reversal 
symmetry and integer angular momentum. It is possible to choose a generic basis, say \(f)k) 
(A; = 1 — > A^), preserving time-reversal symmetry for the matrix representation of H; the 
matrix turns out to be real-symmetric in this basis with its elements H^i = {(j)k\H (pi) . For 
notational simplification, let us denote them by H^ where /i = {kl] s} is a single index which 
can take a value from 1 — > M (M = N{(3N — (3 + 2)/2 the number of independent matrix 
elements). Here (3 is the number of commponents of if^; thus /? = 1 for the real-symmetric 
case. 

Due to presence of complicated interactions in the system, it is technically difficult to 
evaluate some/all elements of the operator matrix in a generic basis. Consequently the 
matrix elements can be determined only within a certain degree of accuracy which, being 



sensitive to local system conditions, varies from element to element. The accuracy fluctuates 
rapidly as system conditions change, with different "time-scale" of fluctuations for each ma- 
trix element. The variation of an element H^ with changing system conditions can therefore 
be mimicked by a particle undergoing Brownian dynamics due to rapidly fluctuating forces 
in addition to an external force (due to existing system conditions). The matrix elements 
of a physical system also have a natural tendency to oppose the cause for their change. The 
dynamics is therefore subjected to a local frictional force too. 

Consider the "particle" H^ in equilibrium under external force V{H^) due to existing 
system conditions. The equation of motion for H^ due to changing system conditions can 
be written as 

d?H„ AH,. 



M _ f^-^-^ti 



Ml S^^^W^Ait,) (1) 

where / is the friction coefficient and A{t^) is a rapidly fluctuating force in "time" t^ 
(a pseudo time only, a measure of the scale for accuracy fluctuations) with following usual 
properties: 

(A(V)A(t^2)...A(t^(2n+l))) = 0, (2) 

{A{t^i)A{t^2)...A{t^i2n))) = Y,^A{t^^A{t^^)){A{t^u)A{t^i))...., (3) 

pairs 

{A{t,{)A{t,,)) = {2/f)6{t,.-t,,). (4) 

where (.) refers to ensemble average, t^j refers to j*^ step in "time"-scale t^ and the sum- 
mation in eq.([n]) extends over all distinct ways in which the 2n indices can be divided into 
n pairs. Further, for a clear exposition of the ideas, we consider here the potential \^(i7^) 
as a function of H^ with a single minima. 

The Langevin equation can now be integrated: let iJ^ be the position of the particle at 
"time" tfj, which changes to position H^ + 6H^ at a later "time" t^ + (5t^ (here t^ chosen 
to be long enough for the effects of initial velocity to become negligible). Due to presence 
of rapidly fluctuating forces, the variation SH^^ in position of the particle will behave like a 
random variables. Using eqs.([T]|l]) and keeping terms only of flrst order in (5t^, one gets 



f{6H,) = -ViH,)St„ f{{5H,f) = ig,/P)St„ (5) 

with g^ = gki = 1 + Sm- Due to random variations in particle position with changing system 
conditions, it is appropriate to consider a "time" -dependent probabihty density p^(if^,t^) 
that the particle will be at the position Hf^ at "time" t^. Assuming a Markovian process 
(that is the independence of future evolution from past states, dependence only on present 
state), one can write 

pAH,, V + K) = / P.(ff„ - SH,; t,)p^AH, - SH,, SH,, St,) iSH, (6) 

where Pcond is the conditional probability that the position of the particle changes from 
Hfj^ — 6H^ to Hfj^ in a time interval 5t^. Expanding both sides of eq.© in a power series of 
6H^ and 6t^ and subsequently using eq.Q we get (in limit 6t ^ 0) 



f 



dp^, _ d 



dt, OH, 



^^ + V{H,) 



Pm (7) 



2/3 dH^ 

Equation([7]) describes the evolution of H^ with respect to "time" -scale t^ which in turn 
depends on the "time"-scale for accuracy-fluctuations (and therefore system conditions) sur- 
rounding H^. For systems where the coupling of any two basis states through the generator 
H is independent of coupling between other states (i.e all matrix elements are independent 
of each other), the fluctuations in accuracy of each matrix element are independent too. 
Each element can therefore be assumed to be subjected to a random force fluctuating at a 
time-scale independent of others, (that is, all t^ independent of each other). This gives us 
M equations, of type (jTj), for the independent evolutions of M elements H^. 

The joint probability distribution p{{H^}; {t^}) of all matrix elements can now be defined 
as 

pi{H,};{t,}) = llp,iH„t,) (8) 



which along with eq.(I7j) leads to the equation for multi-parametric evolution of p: 



dp 



fji^-j: 



d 



OH,. 
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2pdHf, 



+ V{H, 



P 



(9) 
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For a system undergoing evolution as a whole unit, it is natural to seek a common scale, 
say r, at which all its constituents i.e matrix elements vary simultaneously. Let us therefore 
consider the evolution of p with respect to r. Assuming again a Markovian process, we have 



pi{H^}; T + 6r)= f p{{H^} - 6{H^}; r)p,,U{H,} - 6{H^}; 6{H^}; 6r) D6H (10) 

where D6H = Hu'^'^-^m Expanding both sides of eq. lTTOl) in a power series of 6H^ and 6t, 
we get (in limit 6t —>■ 0) 



dp ^^ d 



Or 



dH„ 



dH„ 2 



m,) 



(11) 



As both eq. (ITT]) and eq.([n]) describe the evolution of the probability density of H, they should 
be analogous. A comparison of the equations then gives the conditions 



and 



dp 



6i 



^P'^ 

M ^ 



(12) 



f{6H^) = -V{H,)6t, 
f{{5H,f) = {g,/f3)6T, 



(13) 



The two conditions imply 



St = 6ti = 6t2 



6t 



M- 



(14) 



The above is satisfied if r is defined as r = ^^ '' '" with a„ as arbitrary constants. However 
physical reasoning (based on no preference by random forces to any particular component 
of the system) suggests us to choose Oj equal. 



The solution of eq. ( jTTj) for arbitrary initial condition, say Hq at t = tq can be given as 

p{H, t\Ho, To) oc exp[-{a/2(3)Tr{H - r,Hof] (15) 

with a = (1— r]^)^^ and r] = Q~'y'^^'^°y f . The probability density of if can now be extracted by 
integrating over an ensemble of initial conditions. Although eq. (ITT|) and eq. (IT5l) are derived 
for the case /5 = 1, it is easy to show, following essentially the same steps, their validity for 
the complex Hermitian case {3 = 2 and real-quaternion case /? = 4. 

Note the accuracy scales r^ depend on local system conditions which can vary from sys- 
tem to system. However as eq. lfTT]) and eq. llTSl) indicate, p{H) is insensitive to the details 
of the local system conditions; it depends only on their average behavior described by r 
besides global constraints e.g. V{H) and symmetry conditions; (Note V{H) has no explicit 
dependence on r^). Thus, analogous to their maximum entropy models, the accuracy based 
approach indicates a single parametric dependence of the density p{H) for simple harmonic 
confinement V{H) = H. It further generalizes the formulation to the systems with con- 
ditions giving rise to a generic single-well (single minima) potential. The approach can in 
principle be extended to the multi-well potentials too however it requires a modification of 
the technical details. We intend to pursue these cases in near future. 

It is important to note that the form of eq. OllI) for case V{H) = H is analogous to Dyson 



Brownian model [^, |23l]. The latter deals with the case of a stationary ensemble subjected 
to a random perturbation. However the Brownian dynamics of matrix elements in accuracy 
model is different from Dyson's case; there are two main differences: 

(1) In Dyson's model, the randomness caused due to a perturbation is same for almost 
all matrix elements. In accuracy model, the origin of randomness is the lack of accuracy 
which is sensitive to local conditions. Different matrix elements therefore may be subjected 
to different randomness. 

(2) In Dyson's model, the evolution occurs due to a variation in the perturbation strength 
and is single parametric. In accuracy model, the evolution is brought by the fluctuating 
accuracy due to varying system conditions. As a consequence, we need to consider a multi- 
parametric evolution of probability density (unlike single parametric evolution in Dyson's 
case). However, as eq. flTTl) indicates, the multi-parametric evolution can be reduced to a 
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single parametric evolution. 

III. COMPARISON OF ACCURACY BASED APPROACH AND MAXIMUM EN- 
TROPY APPROACH 

The objective of this section is to indicate the analogy of the results obtained by the 
accuracy model and maximum entropy models of complex systems notwithstanding their 
seemingly different origins. For a clear comparison, we briefly review the maximum entropy 
approach. This approach is based on the representation of a complex system by an ensemble 
of matrices; here the probability density of the matrix elements is formulated by maximizing 
the information entropy under known system-constraints (see [l8| for details). However 
the density in accuracy based approach is obtained as a non- stationary state of a diffusion 
process. This can further be clarified by an example. The accuracy model leads to a Gaussian 
density p{H) if V{H) = H and the initial density is Gaussian too (see eq.( TT5l) . However 
the maximum entropy theory leads to a Gaussian density if the available information about 
matrix elements is limited to their average behavior and variances only: 

p{H, u,b) = l[ p^{H„ v„ h,) = Cexp[- J2i^/2v,)iH, - b,f] (16) 

with C as a normalization constant, v, b as the matrices of variances f^ and mean 6^, 
respectively, and the symbol ^ implying a summation over independent matrix elements 
only. 

The emergence of single parametric formulation in maximum entropy approach can briefly 
be explained as follows. The Gaussian nature of p^ (see eq. (fT6ll ) leads to a relation among 
its derivatives with respect to H^,v^,b^: 

T^p = L^p (17) 



where 

dp dp 



T 



{2/~g^)x^—- - b^ 



dVf, db^ 



L - ^ 



9t^ d ' 

[2/5 SiJ,'"^. 



(18) 
p (19) 



with Xfj, = 1 — g^v^ with g^ = cjki = 2 — 6ki and g^ same as in eq.©. A particular 
combination of the parametric derivatives T = Y2 T^ leads to a diffusion equation Tp = Lp 
(with L = ^^^Lfj,). The single parametric formulation of the diffusion then follows by 
showing T = ^ with Y as the complexity parameter [l9|, l20l |. 

Eq.( TT71) describes the evolution of p when all other matrix elements except iJ^ is held 
fixed. It is therefore equivalent to eq.© with V^Hf^,) = H^ and p^ replaced by p (following 
eq.®. This implies 

The scale t^ can then be expressed in term of the distribution parameters: 

v^ = (1- e^'^'^)/g^ + ci^, 

h^ = e'-/f + C2^. (21) 

with Ci^ and C2^ as constants specific to each f^ and 6^ respectively. The above indicates 
the equivalence of r = Xlu^/^ (^^^ average scale for accuracy fluctuation) to complexity 
parameter in maximum entropy model (an average distribution parameter of the ensemble). 
This further implies that the confinement by a simple harmonic force in the accuracy model is 
equivalent to the maximum entropy modeling of a system with known averages and variances 
of the matrix elements. Similarly a general confining potential V{H) in accuracy model can 
be shown to be equivalent to a maximum entropy ensemble derived under constraints {U{H)) 
=constant where U{H) = J V{H)dH. 

The equivalence of the accuracy approach to maximum entropy approach can be used to 
clarify some of the points related to latter. For example, in maximum entropy approach, a 
particular combination Tp of the parametric derivatives leads to Brownian type diffusion; the 
reason to consider such a combination are not so obvious. However the accuracy approach 
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clearly explains the reason: the combination is required to study the evolution of the system 
as a whole unit. Further, in the maximum entropy approach, the multi-parametric diffusion 
governed by the parameters Yj, J = 12, ...A^ is reduced to a single parametric formulation 
by showing that all Y's except Yi are constants of evolution. However, in accuracy based 
approach, the single parameter existence follows from necessity of the simultaneity of the 
dynamics of various matrix elements. As both approaches represent the same dynamics, this 
reconfirms the lack of any role played by the parameters Y2.., Ym in the diffusion of matrix 
elements. 

IV. DIFFUSION OF EIGENVALUES AND EIGENFUNCTIONS 

The eigenvalue equation of a A^ x A^ Hermitian matrix H is given by HU = U'^E with 
E as the N x N diagonal matrix of eigenvalues, E^n = ^n^mn and U as the N x N eigen- 



vector matrix, unitary in nature: U .U =1 2|]. As described in 211], the statistics of the 



eigenvalues and/or eigenf unctions of H can be obtained from eq. (IT8l) by integrating over the 

211 ] however are valid only for V^ if) 



eigenfunctions, eigenvalues respectively; the results in 
as a simple harmonic force. Here we apply 2nd order standard perturbation theory [27|] to 
derive results for a more general form of V{H)] here again H is taken to be a real-symmetric 
matrix for simplification. 

A. Eigenvalue Statistics 

A small change 6t in the parameter r changes p{H) and its eigenvalue statistics. By 
considering matrix H + 6H in the diagonal representation of matrix H, the change (5e„ in 
the eigenvalues can be given as 

6en = 5H^n + Y. ^^^ + «((^^™")') (22) 
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where Hmn = Gn^mn at value r of the parameter. This further gives, 

N 



f{6en) = {5H^n) + Y. ^^^^^ 



N 



-V{e^)+ Yl 



■^n ^m 



Cn (^m 



6Y 



(23) 



(24) 



Here eq.flMl) has been obtained from eq. fl22l) by using eq. flT3l) . Similarly, upto first order of 
St, 



f{5eje^) = {SH^JHmm) = (2//3)<5„^<5r (25) 

The information about moments of the eigenvalues e„ can now be used to obtain their 



evolution equation. The theory of Brownian motion 



informs us that the joint probability 



distribution P({e„}) for the eigenvalues e^ evolves with increasing r according to Fokker- 
Planck equation, 



Pf 



dP 

Le 



LeP 

N 



E 



n=l 



d 



dCr^ 



_d_ 



+ 



N 

E ■ 



P 



+ men) 



p 



(26) 
(27) 



The above equation describes the evolution of the eigenvalues of a complex system modeled 
by the ensemble p{H) due to changing system conditions. 



As in the case of maximum entropy approach 



19 



22 



2J], the eigenvalue correlations for 



the case V{e) = e can be obtained by using the analogy with Dyson's Brownian ensembles 
[21, l23|]. For a general V{e), the correlations can be analyzed by mapping the eq. (l26l) to 
Calogero-Sutherland Hamiltonian [24]. This can be achieved by using the transformation 



^ = P/\Qn\'^^ in eq.(l26]) reducing it in a form |f = -H^ with Qn = ^„^^n(em - 
e„)'^exp[— /? ^^ f/(e„)] and [/(e) = f deV{e). The 'Hamiltonian' H turns out to be the 
Calogero-Sutherland Hamiltonian in one dimensions 24] 



Y- ^' v:^ /3(2-/9) ^^^,,, 



(28) 



i,j;i<j 
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Similar to the case V{e) = e (see [l9|, l24l for details), the "state" ip or P({e},r|ifo) for a 
generic V{e) can be expressed as a sum over the eigenvalues and eigenf unctions of H. The 
integration of the sum over the initial ensemble Hq would then lead to joint probability dis- 
tribution P{{e},T) and thereby density correlations -R„ for unfolded spectrum (eigenvalues 
rescaled in units of local spectral density). Note the choice of initial eigenvalue distribution 
at To depends on the global system constraints. 

As eg. (1271) and eq. fl28|) indicate, the confining potential V{e) does not affect the short 
range level-correlations. The latter are governed only by the complexity parameter r and 
underlying exact symmetry conditions. However the long range level-correlations are sensi- 
tive to two factors: (i) the complexity parameter (i.e the average accuracy fluctuation scale 
r, or equivalently, the average distribution parameter), (ii) the global system constraints i.e 
details of external force F{H) and the symmetry conditions. Thus the systems subjected 
to effectively similar physical constraints will show analogous long-range correlations (after 
spectral-unfolding) if their complexity parameter are equal. Here the term " effectively similar 
physical constraints" implies the similar symmetry conditions as well as same mathematical 
form of the external potential although it may originate from different physical conditions. 
For example, a harmonic confinement of the matrix elements which also corresponds to their 
Gaussian distribution can be a physical characteristic of many systems related to different 
areas of physics. 

B. Eigenfunction Statistics 

The evolution equation for the probability density of various eigenfunction components 
can similarly be obtained. Here again we consider the case of a real symmetric operator 
for simplification. The eigenvector matrix f/ = O is then orthogonal: 0^.0 = 1 [2]. Us- 
ing standard perturbation theory for Hermitian operators, the second order change in j*^ 
component Ojn of an eigenfunction 0„ due to a small change 6t can be described as 
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^HjnriW^Hm'n ^ 


1 ^ 


- era') 

SHmn 
671 CjnJ 



E \^Hmn\\^Hnn\ ^ _ i_^ ST^ \Ullmn\ /r,Q\ 

As eq.f llSp indicates, the matrix elements of H are uncorrelated. Furthermore, at r, 
-f^mn = Cn^mn (due to H + 5H being considered in diagonal representation oi H) which gives, 
following from eq.([I5]), 5H,nn = -V{Hmn)ST = -V{enS.mn)ST. Thus 6H.mn = V{0)6t = 
for m ^ n and V{0) = 0. The ensemble averaged Ojn then has a non zero contribution only 
from the last term of eq. (!29|) (see eq.Q): 



1 ^ n 

Note for cases where V{Hmn) is nonzero for m y^ n, the first term contributes too. Further 
for cases where V^(0) 7^ or matrix elements are correlated, the other terms may also 
contribute. 

The 2""^ moment of the eigenvector components has a contribution only from the first 
term in eq. fl29l) (up to first order in 6t) 

{60,J0kn)= 2^ (e -e )(E -E ,) ^^"-^^'"' = Z^ (e -e V ^ ^ 

m,m'=pn m=l,m^n 

As the moments for eigenfunction components depend on eigenvalues too, we can first 
write the diffusion equation for the joint probability density Pef,ev{,ei, 62, .., e^; Y) of all the 
components of an eigenfunction and all eigenvalues: 



C-* e/,ei) 



[Lq + LE)Pef,ev (32) 



dr 

where Lq and Le refer to two parts of Fokker-Planck operator corresponding to eigenvalues 
and eigenfunction components. Here Le is given by eq. fl271) and 
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■^o-Y. ' 



dOjn 



[dO, 



^ ((5o,g2) + (5o,„) 



(33) 



A substitution of the moments (eqs.( l30|3T|) ) in eq.( l33l) followed by an integration of eq. (l32|l 
over all eigenvalues except e„ will then lead to the evolution equation for joint probability 



density P„(0„,e„;F); the equation turns out to be same as eq.(18) given in [2l| and can 



further be used to derive various correlation measures for an eigenfunction 2l|. 

The eq.( l29l) can also be used to derive the joint probability distribution of the components 
of different eigenf unctions; again, for the cases with V{0) = 0, the results, e.g. single 



parametric formulation in infinite size limit, turn out to be same as given in 211] 



V. CONCLUSION 

In this paper, we have studied the dynamics of the matrix elements of an Hermitian 
operator of a complex system subjected to a single-well potential. The dynamics is diffusive 
due to random forces originating from accuracy-fluctuations due to varying system condi- 
tions. The information is then applied to explore the statistical behavior of the eigenvalues 
and eigenfunctions. Our analysis suggests a possible classification of complex systems in 
an infinite range of universality classes characterized just by complexity parameter and the 
nature of global physical constraints. The constraints e.g. unitary/ antiunitray symmetries, 
and confining potential on matrix elements seem to divide complex systems in various uni- 
versality classes. Each such class can further be divided into many sub-classes characterized 
by their complexity parameter. Note the "constraint" universality class of a system refers 
to the broad nature of its complexity (the finer details seem to be irrelevant). However its 
sub- universality class depends on the degree of complexity only (measured by complexity 
parameter). This can be explained by following examples. The standard Gaussian orthog- 
onal ensemble (GOE), power law ensemble of real matrices, and, time-reversal Anderson 



ensemble belong to same "constraints" universality class in the above classification [21[ al- 
though their complexity parameters, in general, are not equal (approaching infinity for GOE 
and finite in the other two cases) . However for the system parameters leading to same finite 
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value of the comp 



same statistics 



2l|. 



exity parameter, the Anderson ensemble and power law ensemble show 



The accuracy approach described here is applicable, in its present form, only to the cases 
with independent matrix elements subjected to a single-well potentials. The frequent oc- 
currence of correlated elements or multi-well potentials among complex systems makes their 
analysis desirable too. A generalization to these cases requires a more involved technical 
analysis. However our intuition suggests the possibility of a similar classification for these 
cases too. For example, for the multi-well potentials, the accuracy scales and their fluc- 
tuations are sensitive to local system details and can therefore vary from one branch to 
another. This would lead to a variation of diffusion scales (the average accuracy scale or 
complexity parameter) in different branches. Thus the statistical properties within a single 
branch would belong to a universality class characterized by the local complexity parame- 
ter. However the universality classes in different branches need not be analogous. The above 
suggestion seems to be in accord with already known results for invariant ensembles with 



multi-well potentials |29| . This encourages us to pursue a detailed analysis and extension to 
non-invariant ensembles of such cases in near future. 

For the correlated cases, the accuracy-scales for various elements are no longer indepen- 
dent. However, a recent study of the maximum entropy models of a few correlated cases 
indicates the existence of the universality classes among them too [20]. It is desirable to 
explore the possibility of its generalization to a wider range of such cases. 
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